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1 History and definitions 

The history of Brownian motion was described a number of times. In 1905, 
Einstein established his celebrated formula 

RT 1 



AX* = —- At 

for spherical particles of radius a suspended in a liquid of viscosity /i at tem- 
perature T ; the first member, AX 2 , is the average of the squares of their dis- 
placements during an interval of time At ; R is the constant of perfect gaz, and 
N the Avogadro number. In the following years Jean Perrin made a series of 
experiments leading to a new determination of the Avogadro number, and ob- 
served that the very irregular motion of particles resembled the nowhere diffcr- 
entiable functions of mathematicians. Norbert Wiener introduced what he called 
"the fundamental random function" as a mathematical model for the physical 
Brownian motion. It was called immediately "the Wiener process" , and later on, 
following Paul Levy, "the Brownian motion" . Wiener gave several versions of the 
construction and derived a number of fundamental properties, Levy developed 
the theory to a high point of sophistication, and it is now a mathematical object 
of common use as well as a mine of interesting problems. 

Here is the theory as it appears from the last exposition made by Norbert 
Wiener. The problem is to construct a random process X(t,u>), also denoted by 
X(t) (= X(t, •)), (t the time, u> S fl the probability space) such that 

1) for almost all u> X(t, u) is a continuous function of t 

2) X(t) is a Gaussian process, meaning that the distribution of any n-uple 
(X(ti),X(t 2 ), • • • X{t n )) is Gaussian 

3) this Gaussian process has stationary increments, meaning that the distri- 
bution of X(t) — X(s) depends on t — s only 

4) it satisfies a normalized Einstein equation, that is 



\\X(t)-X(s)\\ 2 2 = \t-s\ 
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where the norm is taken in L 2 (Q). 

Here is such a construction. Let 7i be an infinite-dimensional subspace of L 2 (fi) 
consisting of Gaussian centered variables, and W an isometric linear mapping 
of L 2 (I) (I = K, or 1R + , or [0, 1]) into H. Let xt be the indicator function l[o,t]- 
Then 

X(t) = W( Xt ) 

satisfies all conditions 2) to 4). Moreover, given an orthonormal basis of L 2 {I), 
(u n ), its image by W is a normal sequence (sequence of independent Gaussian 
normalized random variables (£«)), and expanding \t m the form 

Xt = % a n (t)u n (in L 2 (I)) 

results in an expansion of X{t) as a random series of functions : 

X (t) = S a n (t)Cn (in L 2 ([2)) 

or, more explicitly, 

X(t,u) = E a n (t)£ n (w). 

To prove condition 1), it is enough to establish that the series in the second 
member converges uniformly in t for almost all u, and this is done rather easily 
when the u n are classical orthonormal bases. 

By definition, an helix is a curve in a Hilbert space, parametrized by H, 
such that the distance between two points depends only on the distance of the 
parameters : 

\\X(t)-X(8)\\ 2 2 =W-8), 

and ip(-) is called the helix function. A translation of the parameter results in a 
isometric motion of the curve onto itself. It is the abstract model for all Gaussian 
processes with independent increments. 

When ijj(t) = \t\ we say that the curve is a Brownian helix. In contrast with 
the realizations of the Brownian motions (the functions t — ► X(t,u>) when 
lo is fixed), the Brownian helix is a very regular curve. However some basic 
properties of the Brownian motion can be read on the Brownian helix : its 
Hausdorff dimension is 2, its 2-dimcnsional Hausdorff measure is nothing but 
dt, and any three points on the curve are the vertices of a rightanglc triangle 
: the increments starting from a point are orthogonal to the past (therefore, 
independent from the past). 

Simple examples of helices are : 

1) the line (ip(t) = a 2 t 2 ) 

2) the circle (tp(i) — r 2 sin 2 ujt) 

3) the three-dimensional helices (tp(t) — a 2 t 2 + r 2 sin 2 cut) 

4) generalizations of those, with 

ij){t) = a 2 t 2 + J sin 2 art fi(du)) , 

where fi is a positive measure on H + such that the integral is finite. Actually 
this is the general form of an helix function. 



A Brownian quasi-helix in IR 4 built from an automatic sequence 



3 



Except when fi is carried on a finite set, the helix cannot be imbedded in a 
finite dimensional Euclidean space. 

At the end of the 70s, Patrice Assouad developed a theory of Lipschitz em- 
beddings of a metric space into another [2] . He introduced and built quasi-helices 
in Euclidean spaces, meaning that 

\\X(t) - X(s)\\l , 

< a < ^— <b<oo 

ip(t - s) 

for some a and b and all t and s. When tjj(t) = \t\ we call them Brownian quasi- 
helices. Assouad constructed Brownian quasi-helices in Euclidean H™ for n > 3, 
and this gives a new way to prove that the realizations of Brownian motion are 
continuous a.s.. He asked me the question whether a and b can be taken near 
1 when n is large, that is, whether the Brownian helix can be approximated (in 
this sense) by Brownian quasi-helices. I gave a positive answer with an explicit 
construction, and it was published in my paper on Helices and quasi-helices [3] . 



A construction of Brownian quasi-helices by means of 
Walsh matrices 



Let us consider H 2 (n > 1) as a Euclidean space. Let N = 2™. If we want to 
construct a function X : TN — > TR N such that X(0) = and \\X(t) - A(s)|| 2 = 
\t — s\ when \t—s\ < N we have to choose an orthonormal basis uq, u\, . . . ujv-i, 



define 



u N+j 



and write 



x(t)= r 



0<j<t-l 



At this stage there is no restriction on the signs ±, and we may choose + when 
< j < N - 1. If we try to obtain \\X(2t) - A(2s)|| 2 = 2\t - s\, \\X(4t) - 
A" (4s) 1 1 2 = 4|i — s\ etc when \t — s\ < N, whe have more and more conditions on 
the ± and we are led to the following construction. 

We define the Walsh matrix of order N as the N x N matrix obtained as the 

1 1 



,th 



tensor power of the matrix 



1 -1 



that is 



1 1 

1 -1 



1 1 

1 -1 



1 1 

1 -1 



For example, the Walsh matrix of order 4 is 



M = M 2 = 



(11 1 l\ 

1-11-1 
1 +1 -1 -1 

Vi -i-i i j 



(n times) . 
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and the matrix M n+ \ of order 2™ +1 is obtained from M n as 



Af„+i = 



M„ — M, 



) 



The A^ 2 first signs ± are those of the entries of the Walsh matrix, read line by 
line. In order to obtain the following signs, we extend the Walsh matrix by a 
series of vertical translations and change of signs of some lines according to the 
following rule : the first row is nothing but the whole sequence of entries, written 
from line to line and from left to right. 

With this procedure we define X(t) when t is an integer and we can extend 
the construction to all t > 0, then to all real t. It is proved in [3] that we obtain 
a quasi-helix with a and b close to 1 when n is large enough : it is the answer 
to the question of Assouad. 

However, it was not proved that the construction provides a quasi-helix when 
n = 2 (it was remarked that it gives a Peano curve in the plane when n = 1). 
The aim of the present paper is to give a detailed exposition of the case n = 2 
(most of it could be copied for n > 2) and to prove that we obtain a quasi-helix. 
Instead of t £ IR we shall consider only t € IR + and a curve starting from 
(X(0) = 0). We shall investigate the geometric properties of the curve, some of 
them leading to open questions of a combinatorial or arithmetical nature. 

The sequences that we construct are automatic in the sense of [1] 



3 Description of the sequence 



3.1 It is a sequence of +1 and —1 as described before, in case N = 4. We write 
it as a succession of + and — : 

+ + + + +- + - + + -- +--+ + + + + ••• 

The gaps between the blocks of four letters have no meaning, except a help to 
understand the construction. The construction proceeds as follows : given the 
initial word of length 4 J , we divide it into four words A, B, C, D of equal length 
4-?- 1 and write it ABCD ; then 



is the initial word of length 4P +1 . We shall give several equivalent definitions, 
using substitutions, explicit expressions, or generating functions. 



ABC D A (-B) C (-£>) A B (-C)(-D) A{-B){-C) D 
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Beforehand let us write the sequence in a tabular form as in the previous 
section : 



+ - + - 
+ + - - 



+ + + + 
- + - + 



+ + + + 



- + + 



- + - + 

- - + + 
+ - - + 



D 



C 



D 



-B 



a a a t a 2 a 3 
(24 05 



ai6 



a 3 2 



048 



d64 



ai5 



031 



a 47 



063 



079 



3.2 Let me give an explicit expression for a n . Writing 

n = hq + Ani + • • • + A v n v {n v = 1, 2, 3; rij = 0, 1, 2, 3 if j < v) 
the construction shows that 

a n = a-no+irn a m , m = ni + 4n 2 H h A v ~ n v , 

that is 

a n — a no+4ni a ni+4ri2 ' ' ' a ri„_i+4n„ • 

In the second member we find cij s with j < 15. Their value is —1 when j = 
5, 7, 10, 11, 13 and 14 and +1 otherwise. Now let us express a n as a function of 
n written in the 4-adic system of numeration. We obtain the formula 

On - (-l) An 

where A n is the number of 11, 13, 22, 23, 31, 32 in the 4-adic expansion of n. For 
example, ifn=1320011102311122 the significant links are 



1_3_2 1-1-1 2_3_1-1-1 2_2 
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A n is nine and a n = — 1. 



3.3 Let me describe the sequence by means of a substitution rule. 

We start from an alphabet made of eight letters : +a, +b, +c, +d, —a, —6, 
The substitution rule is 



(So) 



The infinite word beginning with +a and invariant under the substitution is 

W = +a + b + c + d + a-b + c-d + a + b-c-d + a-b-c + d-\ 

Replacing a,b,c,d by 1 (or, in a graphic way, in suppressing them), we obtain 
our sequence of ±1 (or ±). 

3.4 Actually there is a simpler substitution rule leading to the same result, 
namely 



+a - 


-> +a + b + 


c + d 


+b 


+a-b + 


c — d 


+c 


+a + b- 


c — d 


+d 


+a-b- 


c + d 


—a 


— a — b — 


c — d 


-b 


-a + b- 


c + d 


—c 


-a-b + 


c + d 


-d 


-a + b + 


c — d 



+a - 


-> +a + 


b 


+b 


+c + 


d 


+c 


+a- 


b 


+d 


+c- 


d 


—a 


— a — 


b 


-b 


— c — 


d 


— c 


-a + 


b 


-d 


-c + 


d 



(Si) 



It can be checked immediately that (Si) (Si) — (So)- 

3.5 The generating function of the sequence (a n ) is 

f(z) = a + aiz + a 2 z 2 + ■ ■ ■ 
It can be defined using partial sums of order 4™ . Let us introduce the matrix 



/I 



M(z) 



1 



z z 
-z z 2 



1 z -z —z' 
\1 -z -z 2 z 3 

and define four sequences of polynomials by the formulas 



[Po\ 






Qo 




1 


Ro 




1 


\T J 




w 



/ Pn+1 
Qn+1 
Rn+l 

V T n+ i J 



= M(z 4 ") 



fP n \ 

Qn 
Rn 

\T n J 
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(n = 0,1,. . .)■ Then 

(Pn\ 
Qn 
Rn 

\Ta) 



M(z 4 " )M{z 4 )---M(z 4 )M{z) 



1 
1 

W 



When \z\ = 1, we have M{z)M(z) — 41, therefore the matrix ^M(z) is unitary, 
hence 



\P n \ 2 + \Qn\ 2 + \R n \ 2 + \T n \ 2 = 4(| J P„_ 1 | 2 + \Q n _^ + |i? n _r| 2 + \T n ^) 

= 4"(1 + 1 + 1 + 1) = 4™ +1 

We obtain the generating function as 

f(z) = lim P„(z) 

n — *oo 

We can write the generating function in a more interesting form : 
f(z) = / (z 4 ) + z/i(z 4 ) + z 2 / 2 (z 4 ) + z 3 / 3 (^ 4 ) , 

where the coefficients of the power series /c/1,/2,/3 are the columns of the 
table in 3.0. In order to obtain these coefficients, we can start from W in 3.2 
and replace a, b, c, d by 1, 1, 1, 1 (for / ), 1,-1, 1, -1 (for /1), 1, 1, -1, -1 (for / 2 ) 
and 1, -1, -1, 1 (for / 3 ). Then / = /. Writing 



F(z) = 



h{z) 



the functional equation of the generating functions of the columns is 

F(z) = M{z)F{z A ) . 



4 Description of the curve 

4.1 Let wo, Mi) u 2, U3 be an orthonormal basis of the Euclidean space H 4 , and 
define v,j +4 = Uj (j = 0, 1, . . .). The partial sums of the series 

a u + a i u 1 + a 2 u 2 H 

(that can be obtained from W in 3.2 by replacing a, b, c, d by u , ui,u 2 , u 3 ) will 
be denoted by S(n). Then 

S(n) = a u a + aiui H h a„_iu„_i e 2Z 4 . 

It is easy to check on the table in 3.0 that 

S(Wn) = 4S(n) (n e IN) . 
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Moreover it is not difficult to see (we shall be more specific later) that 

\\S(n)- S(m)\\ 2 <b\n-m\ 

for some b < oo and all n and m. This allows, first, to define S(t) when t is a 
binary number via a formula S(16"i) = l v S(t), then to check that 

5(16t) = 4S(t) 
||5(t)-5( s )|| 2 <6|t- s | 

for such numbers, then to extend S(-) by continuity on M + , and check the above 
formulas for all t > and s > 0. 

The curve we consider is the image of H + by S(-). 

Clearly (changing t into 16t) the curve is invariant under an homothety of 
center and ratio 4. Our main purpose is to prove that it is a Brownian quasi- 
helix. We shall point out some geometric properties first. 

4.2 The matrix M transforms Uo, u\,u 2 , M3 into u' , u' 2 , u' 3 : 

(u' ,u'i,u' 2 ,u' 3 ) = M(u Q ,ui,u 2 ,u 3 ) 

and the partial sums of order n of the series Saju'j are the partial sums of order 
An of the series EajUj. Therefore the equation 

S(4i) = M Sit) 

holds true when t — n S IN and by extension for all t € IR + . 

It is easy to check that the eigenvalues of M are 2 and —2, and that 



M 



Let 

w = \{u + Mi + u 2 - Us) , Vi = ^-(uq + u 3 ) 

V2 = \{U() - Ml - M 2 - M 3 ) , V 3 = ^{Ul - M 2 ) . 

They constitute an orthonormal system. The vectors i>o and V\ generate a plane, 
P, which is the eigenspace of the eigenvalue 2, and v 2 and M3 a plane, Q, corre- 
sponding to the eigenvalue —2. Expressed via the orthonormal basis vo, v\, V2, V3, 
the operator M takes the form 

^10 \ 
10 
0-10 
\0 -1 J 







( 2 2 ^ 




( 1 


°\ 




-2 


\ 


1 




2 


M 


-1 


1 




2 


-2 


1 




2 


-1 


-1 




2 


2 






\-22) 






/ 


V 


2 


/ 



M' = U M U- 1 = 2 



where J7 is the unitary matrix carrying uq, u\, u 2 , M3 onto vq, Vi, v 2 , V3. It means 
that the transformation S(t) — > 5'(4t) is the product of a homothety of centre 
and ratio 2 and an orthogonal symmetry with respect to the plane Q. 
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4.3 Clearly M 2 — AI (I being the identity matrix), In turn, M is the square 
of another simple matrix, as it can be guessed from 3.2 and 3.3. Let us dehnc 



T 



(l 1 \ 
10-10 
10 1 

\0 1 -I J 



Then M = T 2 . 

The eigenvalues of T are y/2, — y/2, iV% and —i^/2. The vectors 



w 



U>1 



are eigenvectors corresponding to \/2 and — \/2- Dchning w 2 and w 3 in such a 
way that wo,wi,W2, W3 is a direct orthonormal basis, and W being the unitary 
matrix carrying uq, u\, 112, U3 onto wo, w\, W2, W3, we can write 



wtw- 1 = T 1 = V2 



(I \ 
0-100 
-1 

\0 1 j 



It means that T' is decomposed into : 

1) an homothety of centre and ratio \/~2 

2) a rotation of | of the orthogonal projection on Q 

3) a symmetry with respect to w\ of the orthogonal projection on P. 

In the same way as we obtained the equation S{At) — MS(t), we now have 



S(2t) — T S(t) 

and we have just given the interpretation of the transformation S(t) 
as a product of simple transformations. 



S(2t) 



4.4 We have investigated the properties of the transformations S(t) — > S(16t), 
S(t) — ► S(At), S(t) — ► S(2t) as products of homotheties and isometries. Now 
we shall look at the effect of a translation of t by an integer. We are interested 
in differences S(t) — S(s). 

Let us begin with integers m < n < 16 k . Let us divide the series aou + 
aiui + ■ ■ ■ into consecutive blocks of length 16 fc , so that the series reads 

+A + B + C + D + A-B + C-D + --- 

If j = jo + 4ji (j = 0, 1,2,3, ji G IN), the j-th term is of type A,B,C,D 
according to the value of jo and its sign is aj . Therefore 



S(n + j ■ 16 fe ) - S(m + j ■ 16 fe ) = aj (S(n + j 16 fe ) - S(m + j 16 fc )) . 
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If < s < f < 1, we can approximate s and t by m 16 k and n 16 k and we 
obtain 

S(t + j) - S(s + j) = aj (S(t + J0 ) - S(s + j Q )) 

(jo = 1, 2, 3, 4, j = j modulo 4). 

This expresses that all arcs Aj = S(\j, are isometric (actually translates 

or symmetric according to the value of aj) of one of the arcs Ao,Ai,A2,Aa 
(according to the value of jo). Using 4.2, this holds true when we replace the 
^•s by A) = S{\j2 v , (j + l)2 u \) whatever v e ZZ. 

5 It is a Brownian helix 



5.1 What we have to prove is that, writing 

, ||S(t)-S(s)|| ||5(t) -5( 8 )|| 
a = ml ; — < sup . — = o , 

0<s<t v\t-s\ 0<s <* VI* ~ s l 

we have 

a > , 6 < oo . 

We can write as well 



; ™ ^/|n — m\ m<n y/\n — m\ 



5.2 The easy part is b < oo. 

Let us first assume [m,n] = [j2 fc , (j + l)2 fc ]. Then, according to 4.3, 

\\S(n)-S(m)\\ = 2 k /\ 

In the general case, let us decompose [m, n] into such intervals in a minimal way, 
so that there are at most two intervals of the same length in the decomposition. 
If the largest length is 2 k , we obtain 

\\S(n) - S(m)\\ < 2(2 k ' 2 + 2 {k -^' 2 + •••)< 2(1 + V2)2 fc/2 

therefore 

\\S(n) - S(m)\\ <2(1 + V2)| 

This gives 

b < 2(1 + V2) 

5.3 To prove a > is more tricky. 
We shall use two lemmas. 

Lemma 1. There exists a > such that 

\\S(n + h)-S(n)\\ <l-a 

for alln € IN and he [—5, 3] fwif/i n + h > 0). 
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Lemma 2. There exists an integer A such that 

||S(n)-S(ro)||>2 
for all integers m and n such that n — m > A. 

Assuming that this is correct, the result is at hand : given t and s such that 
t — s > A + 1, we can write s = m + h and t = n + h! with h, h' e \—\, \] and 
m — n > A, therefore 

\\S(t)-S(s)\\>2a. 
Whenever (A + l)2 fe <t-s<(A + l)2 fe+1 (k E ZZ), we have 

2a 



\\S(t)- > 2a2 fe / 2 > 



y/2(A+l) 



,1/2 



therefore 



a > 



2a 



y/2(A+l) 



5.4 Proof of Lemma 1 . 

From now on it may be useful to represent S(n) on the table of 3.0, and 
also the differences S(n) — S(m), as figures consisting of consecutive lines plus 
or minus part of a line above and below, in such a way that each column in the 
figure has a sum equal to the corresponding coordinate of S(n) or S(n) — S(m). 




5(17) 




5(32) - 5(25) 



Let us consider 



||5(16n + m)-5(16n)|| 2 

for n € IN and m — 0, ±1,±2, ±3, . . .,±8. It is sufficient to consider the four 
cases n = 0,1,2,3, and to look at the figures (depending on m) in each case. 
The result is 

||5(16n + m) - 5(16n)|| 2 < 8 (n odd) 
||5(16n + m) - 5(16n)|| 2 < 9 (n even) 

with equality only when m is odd (as for 5(32) — 5(25)). Therefore, going one 
step further, 

||5(l6n + m+ - 5(16n)|| < V8+ \V9 (n odd) 
||5(l6n + m+ - 5(16n)|| < V9+ ±V8 (n even 
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when p = 0, ±1, ±2, . . . , ±8. Proceeding that way we finally obtain 
||5(16(n + t)) - 5(16n)|| < 4(1 - a) ( - I < i < I) 

where 

4(l-a) = (V9+JV8)(l + ^ + I ^ + -) 
as the second member is less than 4 and that proves Lemma 1. 

5.5 Proof of Lemma 2 

Here again we look at the table. We can compute | \S(n) — S(m) \ | when n — m 
is in a given interval, and moreover give the couples (to, n) for each the infimum 
is attained, and the expression of 5(n) — S(m) (that is, the coordinates with 
respect to uo, ui,U2,u 3 ). 

1) 4 < n — m < 16. It suffices to consider the first 16 lines of the table (fig. 1), 
since adding to m and n a multiple of 64 does not change \ \S(n) — 5(to)||. We 
obtain 

inf \\S(n) -5(to)|| 2 = 2 
realized for (5,11), (23,29), (35,41) and (53,59) : 



5(11) - 5(5) = m - u 4 , 5(29) - 5(23) = u 2 - u 3 
5(41) - 5(35) = -u 2 + u 3 , 5(59) - 5(53) = -u x + u 4 



5(11) -5(5) 

5(29) - 5(23) 
5(41) - 5(35) 

5(59) - 5(53) 



- + + - 



t 6 
t 7 
£ 8 
I 9 
t 10 
t 11 
5(42) - 5(22) 



fig. 2 



fig. 1 



2) 16 < n — m < 64. It suffices now to consider the first 256 terms (the 
first 64 lines of the table). The idea in order to pick the infimum is to start from 
5(4 x 11) — 5(4 x 5) and the analogues, and modify the figure in order to diminish 
\\S(n)-S(m)\\ 2 (fig. 2). As a first example, 5(44) -5(20) = 2u 2 + 2u 3 (obtained 
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by replacing u\ + u 2 + u 3 + w 4 and u 4 by u\ — u 2 — u 3 + w 4 in the expression of 
5(11) — 5(5)), and the modification provides £(42) — 5(22) = u\ + u 2 + u 3 — u 4 . 
The result is 

inf||5(n)-5(™)|| 2 = 4 
realized for (22,42) and (214,234), with 

5(42) - 5(22) = mi + u 2 - u 3 - w 4 , 
5(234) - 5(214) = -u x -u 2 -u 3 + u A . 

This proves Lemma 2 with A = 16. 

Actually the proof can be given in a more concentrated form. It is enough to 
show that \\S(n) — S(m)\\ 2 < 3 is impossible when n — m > 16. Let us assume ab 
absurdo that \\S(n) - S(m)\\ 2 < 3. Let us add or remove the minimal number of 
terms in order to transform 5(n) — S(m) into a difference of the form 5(4n') — 
5(4m') (that is, to transform the figure S(n) — S(m) into a rectangle). In general, 
this minimal number is < 4 and has the same parity as ||(5(n) — 5(m)|| 2 ; here 
it is < 3 and the resulting 5(4n') — 5(4m') has its squared norm < 2, therefore 
||5(n') — 5(to')|| 2 < 3 and the process goes on until we reach S(n x ) — S(m x ) 
with n x < 64. Then we know the possible pairs (m x , n x ), namely (5,11), (23,29), 
(35,41) and (53,59), and the reverse process never gives a squared norm < 3. 

5.6 Remarks and questions 

The estimates we gave for b and a are quite rough. We can ask for better es- 
timates and conjectures. The actual problem, of a combinatorial or arithmetical 
nature, is to compute these numbers exactly. 

We were interested in estimating b from above and a from below. Examples 
provide estimates in the opposite direction : 

b > mm = * > 1.21 

— VTr vvt — 
a < LM^CT = 2 = 1 < 0.45 

It seems not impossible that the estimate for a is precise, that is a = -^=. 
A careful investigation of the table would confirm or disprove this conjecture. It 
would lead also to a better estimate for b. 



6 Projections of the curve 

6.1 The direction of u$ is special : all first coordinates of the 5(n) are > 0. 
That means that the partial sums 5o(n) of the original series described in 3.0 
are positive. 

A simple way to see it is to use Lemma 1 (of 5.2). Since 5o(n) > 1 for 
n = 1,2,3,4,5,6,7,8, we have 5o(t) > for \ < t < 8, therefore (changing t 
into 16 fe i), S (t) > for all t > 0. 
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6.2 We are mainly interested in the three-dimensional projections of the curve. 
It seems likely that all parallel projections of the curve on a three-dimensional 
subspace of H 4 have an infinity of double points. The question can be formulated 
in the equivalent forms : 

1) is every direction in IR 4 the direction of some S(t) — S(s) ? 

2) are the S( " ) ~ S(m) (n > m) dense on the sphere S 3 ? 

6.3 Let us project the curve from on the sphere S 3 , that is consider 



5W = Pii ( * >0) - 

We obtain a closed C, image of any interval [a, 16a] by S(-). 

C is invariant under the isometries of H 4 defined by and -^T (see 4.1 
and 4.2). The first takes the form 

/l \ 
10 
0-10 
\0 -1/ 

with respect to the orthonormal basis (vo, fi, «2, ^3) and the second 

/l \ 
0-100 
-1 
\o 1 / 

with respect to the orthonormal basis (wo, w\, W2, W3). The vectors vo and vi 
(wo and w\ as well) generate a plane P such that the mapping S(t) — ► S(4t) 
is an orthogonal symmetry with respect to P. For the projection of C on P, the 
change of t into 2t means a symmetry with respect to the line generated by w . 

C has a double point at t = \, if = f : T(|) =T{\). 

In order to prove it, we expand t and t' in base 4 (we underline the expansion) 




V2 



r 



t = 0.1111 • • • if = 1.1111 • 



A Brownian quasi-helix in IR 4 built from an automatic sequence 
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Using base 4 again, we easily obtain the figures and the values of 5(1), 5(11), 
5(111) and so on : 5(1) = u = (1 0) 

5(11) = 5(10) + (5(11) - 5(10)) = (1 1 1 1) + (1 0) 

5(111) = 5(100) + (5(110) - 5(100)) + 5(111) - (5(110) 

= (4 0) + (1 1 1 1) - (1 ) 

5(1111) = 5(1000) + (5(1100) - 5(1000)) + (5(1110) - 5(1100)) 

+(5(1111) - 5(1110)) 

= (4 4 4 4) + (4 0) - (1 1 1 1) + (1 0) 

5( 11111 ) = (16 0) + (4 4 4 4) - (4 0) + (1 1 1 1) - (1 0) 

5( 111111 ) = (16 16 16 16) + (16 0) - (4 4 4 4) + (4 0) 

-(1 1 1 1) + (1 0) 

The ratio between two consecutive vectors tend to 2 (meaning that the ratios 
of coordinates tend to 2), hence 

S(t') = T(t) (t= 1 -). 
By isometry we also have 

T(2i') = T{2t) . 

These double points are contained in the plane P, and they are symmetric 
with respect to the line generated by wq- 

I believe, but did not prove, that these arc the only multiple points of the 
curve C. In that case, C is a Brownian quasi-helix (actually, a Brownian quasi- 
circle) on some 4-covering of the sphere 5 3 . 

6.4 One can see the curve C in two other ways 

First, taking into account that the first coordinate So(t) is always positive, 
we can consider 

nit) M. 

m S (t) 

and the curve C described by IZ(-), projection of the original curve with a source 
at and a screen at the hyperplane x = 1. 

Symmetries and double points can be studied on this model as well. 

Secondly, we obtain a projective model of C, say, C", on choosing four points 
A , Ai, A 2 , A 3 in K 4 , defining A/+4 = A, (j = 0, 1, . . .), starting with a point 
M = A and defining the sequence of points 

M n+1 = ■ ((a + oi H a„_i)M„ + a n A n ) . 

ao + ai H a n 



Some real figures would help. If a reader is willing to draw figures of the 
above curves, I'll appreciate to see them. 
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